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ABSTRACT 

Using a dark matter only Constrained Local UniversE Simulation (CLUES) we examine the 
existence of subhaloes that change their affiliation from one of the two prominent hosts in 
the Local Group (i.e. the Milky Way and the Andromeda galaxy) to the other, and call these 
objects "renegade subhaloes". In light of recent claims that the two Magellanic Clouds (MCs) 
may have originated from another region (or even the outskirts) of the Local Group or that they 
have been spawned by a major merger in the past of the Andromeda galaxy, we investigate 
the nature of such events. However, we cannot confirm that renegade subhaloes enter as deep 
into the potential well of their present host nor that they share the most simplest properties 
with the MCs, namely mass and relative velocity. Our simulation rather suggests that these 
renegade subhaloes appear to be flying past one host before being pulled into the other. A 
merger is not required to trigger such an event, it is rather the distinct environment of our 
simulated Local Group facilitating such behavior. Since just a small fraction of the full z = 
subhalo population are renegades, our study indicates that it will be intrinsically difficult to 
distinguish them despite clear differences in their velocity, radial distribution, shape and spin 
parameter distributions. 

Key words: cosmology: theory - cosmology: dark matter - methods: A-body simulations - 
galaxies: Local Group - galaxies: Magellanic Clouds 
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1 INTRODUCTION 

In the concordance cosmology, structure forms in a hierarchical, 
"bottom-up" fashion that leads to the accretion of small substruc- 
tures by large dark matter haloes. The current paradigm holds that 
substructures orbit within their parent haloes until tidal stripping 
rips them apart as they sink to the centre of the potential by dynam- 
ical friction. 

The Local Group o f galaxies is an ideal testbed for th is kind 
of near-field cosmology JFreeman & Bland-Hawthorn||2002l) as the 
observational data is better than for any distant galaxy. But as al- 
ready prompted before, how certain can we be that the Milky Way 
(MW) or the Andromeda galaxy (M 31) are in fact typical galax - 
ies of their mass or luminosity (cf. iForero-Romero et alJuOl lh ? 
These two galaxies are the dominant members of the Local Group 
and form a two-body system on a collision course i n approximately 
2-3 Gyrs dCox & Loebll2008l : lHoffman et alj|2007l) . Their proxim- 
ity at the present time indicates that they may already have started 
to influence each other or their respective satellite populations. Per- 
haps thjsjnfluencejs related to the (in)famous "disk of satellites" 
(e.g. lMetz e t al. 2008); there is also the question of how frequently 
the largest companions of the MW - the Small and the Large Mag- 
ellanic Clouds (SMC and LMC) - are found. A recent study by 



iLiu et alj 1 120 101) indicates that only of order 3% of Milky Way type 
galaxies in the SDSS Data Release 7 host two satellites with lumi- 
nosities simi lar to the Magellanic C louds (MC). This result is also 
supported bv llames & Ivorvl ( 120101) who investigated 143 luminous 
spirals in Ha finding that the MW is an unusual galaxy both for the 
luminosity and the proximity of its two brightest satellites. There- 
fore, both these satellites are often considered outliers with respect 
to the full system of satellites or biting the MW and their forma- 
tion scenario is a matte r of debate dBesla et alj2007l : |Peeblesl2009i 
iKallivavalil et al.ll2009l : lMetz et alj2009h . 

There are several investigations that suggest that the LMC and 
SMC are on their fi rst infall as sugge s ted by recent proper mo- 
tions measurements dBesla et alj 120071: iTollerud et alJuOl lh - as 
opposed the to recent claims by ISales et alj ( 1201 II) that the LMC 
is not necessarily on its first approach to the MW based upon the 
Aquarius simulation dSpringel et al.ll2008l) : if they are only loosely 
(if at all) bound satellites of t he MW then where did they come 
from? lYang & Hammen ( 120101) deal with the possibility that both 
the Magellanic Clouds have been expelled from M3 1 due to a ma- 
jor merger event. The scenario envisioned by them is that they are 
ejected tidal dwarf gala xies from a previous major merger occur- 
ring at the M31 location JHammer et al.l2010f) . But there is also the 
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notion that they may simply be a ccreted objects that formed in the 
outer reaches of the Local Group jvan den Berghl201(j| : lBusha et al.l 
120101) . While these studies are still speculative, they nevertheless 
show that M3 1 may play a vital role in shaping the orbits of the 
Magellanic Clouds (see e.g. iKallivavalil et al.ll20(M) . Neither the 
MW nor M3 1 may be understood in isolation but only in the context 
of the (unique) environment of the Local Group and its formation. 
In this Letter we use a Constrained Local UniversE Simulation 
(http://www.cIues-project.org) to search for events where subhaloes 
change their host halo affiliation. These simulations are well suited 
for this idea as they directly model the Local Group (consisting of 
the two-body system MW and M3 1 ) within the correct environment 
and in a cosmological framework of the concordance model. We 
investigate the likelihood that subhaloes were under the influence 
of one of the two hosts at some previous time while at z — are 
within the virial radius of the other host. We colloquially call them 
"renegade" or "disloyal" subhaloes. 
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2 THE SIMULATIONS 

Constrained Simulation of the Local Group We choose to run 
a dark matter only simulation using stan dard ACDM in itial condi- 
tions, that assume a WMAP5 cosmology (Komatsu et al.l2009f) . i.e. 
Q m = 0.233, il b = 0.046, Cl A = 0.721. We use a normalization 
of us = 0.817 and an = 0.96 slope of the power sp ectrum. We 
use the treePM-SPH code GAD GET 2 JSpringelll2005l) to simulate 
the evolution of a cosmological box with side length of Lbox = 
64/i~ Mpc. Within this box we identified (in a lower-resolution 
run utilizing 1024' ! particles) the position of a m odel local group 
that c losely resembles the real Local Group (cf. iLibeskind et al.l 
I2OIO1) . This Local Group has then been re-sampled with 64 times 
higher mass resolution in a region of 2/i _1 Mpc about its centre 
giving a nominal resolution equivalent to 4096 3 particles giving a 
mass resolution of rr tpM = 2.1 x 10 5 h~ 1 M p ) . For more details we 
refer to the reader to lGottlober et al.l J2010h . 



The (Sub-)Halo Finding and Trackin g We use the 

MPI+ OpenMP hybrid halo finder AHF JKnollmann & Knebe] 
120090 to identify haloes and subhaloes in our simulation which is 
the successor of the MHF halo finder by I Gill et all J2004I) . We only 
consider subhaloes where the number of gravitationally bound 
particles is larger than 20. 

We build merger trees by cross-correlating haloes in consecu- 
tive simulation outputs. The direct progenitor at the previous red- 
shift is the object that shares the most particles with the present 
halo and is c losest to it in mass. Fo r more details we poin t to the 
reader to e.g. ILibeskind et alj J20T0I) or lKnebe etal] J20ld) . 



Seeking Renegade Subhaloes The prime focus of this Letter is 
to examine the set of subhaloes that change their host affiliation. 
These subhaloes are identified as being within a fixed distance of 
one host at some early time and within the same fixed distance of 
the other host at z = 0. In practice we identify renegade subhaloes 
within D = i?5° r st of the two main members of the Local Group 
(MW and M3 1) at redshift z — 0. Each subhalo is then traced back 
in time using the merger tree. At each snapshot we compare its 
distance to the other host and note when this distance falls below 
D{z) = Rvi^iz), i.e. when it is accreted, recording it as a "host 
change". 



Figure 1. Projection of the nights paths in comoving coordinates of both 
hosts MW (thick dashed) and M31 (thick solid) as well as those satellites 
that changed loyalty. The crosses mark the subhaloes' staring points at high 
redshift whereas the end points of 2 = are unmarked. Please note that we 
only show the trajectories for subhaloes with more than 1000 particles. 



3 RESULTS 

We begin this section by describing how disloyalty in subhaloes 
may arise. In Fig[T]we show renegade subhaloes of our simulated 
local group. The two thick lines represent the flight paths of the 
"MW" (dashed) and "M31" (solid) from z = 5 until today, from 
now on simply referred to as MW and M3 1 : they are currently ap- 
proaching each other with 113 km/s, having a (relative) transversal 
velocity of 23 km/s. In blue and red we show the trajectories of 
renegade subhaloes of the simulated MW and M3 1 , respectively. 
Please note that this figures serves merely an illustrative purpose 
and hence has been limited to just the most massive subhaloes, 
those with more than 1000 particles. Including lower mass sub- 
haloes has no qualitative effect. We actually observe the same phe- 
nomenon reported and quantified for the CLUES simulation within 
a WMAP3 framework: (renegade) subhaloes fall in to their respec- 
tive (first) hosts from preferred directions, i.e. they are residing in- 
side the filamentary structure throug h which the two host s travel 
towards the (simulated) Virgo cluster JLibeskind et alj201 ll ). How- 
ever, while this filament is broken into two strands in the WMAP3 
simulation, a visual impression (not presented here) of the WMAP5 
simulation reveals that the MW and M31 are occupying the same 
filament (Libeskind et al., in preparation). This explains the possi- 
bility to exchange subhaloes whereas we cannot confirm their exis- 
tence in the WMAP3 simulation (see below). 

We find that the situation is slightly asymmetric with respect 
to the two hosts: there are 128 renegades in total: 72 of the MW 
and just 57 of M31 when considering only well resolved subhaloes 
(i.e. those with more than 20 particles) found by our halo finder 
AHF. Note that we fail to find multiple host changes: if a subhalo 
changes its host it happens just once (or never) in our particular 
simulation. The number of renegade subhaloes (and their fractions 
with respects to the total number of subhaloes) is summarized in 
Table [1] 

We apply our algorithm at higher redshift when the two haloes 
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Table 1. The number of renegade satellites (with N p , SSjt ^ 20) as a func- 
tion of redshift z and distance £>mw-M31 ( as measured in h~~ 1 kpc), re- 
spectively. The last line summarized the percentage p of renegade subhaloes 
with respects to the total number of subhaloes. 
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were further apart; note that the two hosts monotonically approach 
each other since z ~ 5. We find that the number of renegade sub- 
haloes is a function of host distance - the closer the two haloes are 
to each other, the more renegades exist. We thus conclude that the 
number of renegades is an indicator of the continuous approach of 
the MW and M3 1 and foretells their imminent merger, respectively 
JHoffman et alj|2007l) . We repeated our analysis in the WMAP3 
CLUES simulation (e.g. lLibeskind et alj|20ict iKnebe et alj|201lh 
where the 2 = separation of the two hosts is ~ 700/i _1 kpc. In 
this case we find no renegades. We ascribe the decreased number to 
the lower erg -normalisation of that particular model and the smaller 
masses of the two respective hosts, respectively. Further, whereas 
the two hosts appear to be flying within one common filament in 
the WMAP5 simulation only gradually moving towards each other 
over time, the orbital approach in the WMAP3 simulation is more 
rapid close to redshift z = due to the embedding of each host in 
a separate filament (cf. Libeskind et al., in preparation). 

While Fig[T] shows the trajectories of the renegade subhaloes 
it does not provide information about the first (host) infall, exit and 
second (host) infall times. To gain insight into this issue we present 
in Fig[2]the corresponding redshift distributions for the combined 
sample of the MW and M31: a peaked dis tribution can be in - 
dicative of g roup infa ll as reported before by Li & Helrnil 020080 : 
lAngulo et all < J2009t) ; [Klimentowski et alj tepid) . Here we find that 
even the combined sample shows a peaked distribution meaning 
that infall and exit happened for both populations of renegades (i.e. 
MW and M3 1 renegades) at approximately the same time. We find 
that (most of) the renegades fell in at approximately z (a 0.55, 
left it shortly afterwards (z ss 0.35) to enter the other host close to 
redshift z = 0, This confirms the picture already dr awn from Fig[J 
reneg ade subhaloes come from the same region (cf. lLibeskind et al.l 
1201 If) , fly past one host, and are (gravitationally) pulled over to the 
other host, ending up within the area of trade of the latter at redshift 
z = 0. It is importa nt to note that despite claims by, for example, 
iFouquet et alj J201ll) that a major merger is required to trigger rene- 
gade satellites, we find them without such events. In our simulation 
the constantly decreasing distance and proximity, respectively, be- 
tween the two hosts and the existence of backsplash galaxies (i.e. 
galaxies that enter and leave the virial radius of a host halo; cf. 
iGill et aljUooH ; IKnebe et al1l201ll) is sufficient to explain the ap- 
pearance of renegade subhaloes. 

The next question we address is that of mass distribution of 
the renegade subhaloes. The corresponding plot is to be found in 
Fig[3]where we show the mass spectrum of the disloyal subhaloes 
at redshift z = (solid line) in comparison to the full subhalo 
mass function (dashed line), as an average of the MW and M31. 
We additionally show the distribution of renegade halo masses at 
the time of the initial infall into their first host as dotted line. 

We observe that both mass functions at redshift z — follows 
a power-law. Or in other words, the renegade function is not peaked 
at a particular mass. However, the masses of the renegades are sys- 
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Figure 2. The distribution of first infall, exit and second infall redshifts for 
all renegade subhaloes. 
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Figure 3. The mass function of renegade subhaloes at redshift z = (solid 
line) and at their respective infall time into the first host (dotted line) com- 
bined for M31 and the MW. We also show the (combined) mass function 
for all present-day subhaloes (dashed line). 



tematically lower than those of their loyal companions. Further, it 
is also apparent the renegades are less massive than the LMC yet 
compatible with the SMC. Our simulations therefore provide evi- 
dence for the existence of renegades when the environmental set- 
ting of the hosts is just about right, but the link to the Magellanic 
Clouds is not immediately apparant. That being said, our inabil- 
ity to produce MC-like renegade subhaloes may just be due the 
difficulty in forming such massive satellites in the fi rst place (e.g. 
iBovlan-Kolchin et al.|[2oTol J201 ll : JSales et alj|201lh . Figj3] further 
confirms that the renegades suffered mass loss as their original 
masses at first infall were around a factor of two larger. 

We additionally query whether the unique past of renegade 
subhaloes will leave a noticeable imprint upon their (internal) prop- 
erties other than their mass. To this extent we calculated several 
quantities of interest including their relative velocity (within the 
rest frame of their respective host), their distance to the final host 
at redshift 2 = 0, their sphericity s = c/a based upon the largest 
and smallest eige nvalue (a and c, re s pectiv ely) of the moment of 
inertia tensor (cf. IKnebe et alJuOlfJ. 120081) . their spin parameter 
A = \L\I\/2MVR JBullock et alj|200ll) . and their concentration 
as defined by x = r su b/r2 (where r su b is the subhalo edge and 
T2 the peak position of p(r)r 2 with p(r) being the density profile). 
For this particular test and any possible comparison to the Mag- 
ellanic Clouds we restricted all subhaloes to contain at least 500 
particles leaving us with 18 renegades and 206 subhaloes in total. 
The solid lines in Fig|4] refer to the distribution of the renegade 
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Figure 4. Distribution of relative velocities (upper left), distance (upper right), sphericity (lower left), and spin parameter (lower right) for renegade subhaloes 
(solid) in comparison to all (dashed) and loyal (dotted) subhaloes at redshift 2 = combined for the MW and M3 1 . 



subhaloes alone, while the dashed line represents the total subhalo 
population (including renegades); the dotted line represents just the 
non-renegade objects. Each curve is normalized by the respective 
total number of subhaloes. 

From Fig|4]we see that renegades posses substantially larger 
2 = velocities. However, hardly any of the renegades show a ve- 
locity as large as the one observed for the Ma g ellanic Clouds (e 
340-380 km/sec, see iKallivavalil et alj|2006l ; IPiatek et alj|200: 
likely because they do not enter as deep into the potential as their 
observational counterparts. Furthermore the most massive rene- 
gades are not the ones with the largest relative velocities. They also 
do not enter as deep into the potential as their loyal companions. 
However, we need to acknowledge that our simulation does pro- 
duce non-zero objects with velocities comparable to the MC's even 
though their number is rather low. When examining the sphericity 
and spin of each renegade, we find that the they appear rounder and 
with a higher spin parameter than the general population. We ex- 
amined the distributions of concentrations and found them to be in- 
distinguishable from the full subhalo population (not shown here). 
Despite the apparent differences seen in Fig|4]we nevertheless con- 
jecture that they are unable to leave an discernible imprint on the 
respective distributions due to their low numbers. 

We finally examine the spatial distribution of our renegade 
subhalos by showing, in Fig[5] the distribution of cos 8 = R sa t • 
Dmw-m3i where R sa t is the (normalized) position vector of a 
subhalo in the frame of its redshift z = host and Dmw-m31 
is the (normalized) vector connecting the MW to M31. We can 
clearly see that this angle is isotropically distributed for the the 
full subhalo population of the two hosts, while, there is a clear 
preference for renegade satellites to be found in the direction of 
the other host. Please note that the near uniform distribution about 
Dmw-m3i found for the loyal satellites here does not contra- 
dict the generally accepted anisotropic distr i bution of subhaloes 
(e.g. iKnebe et alJl2004HLibeskind et al. II 20051 : IZentner et ai]|2005t 
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Figure 5. Distribution of the (cosine of the) angle between the position 
vector of a subhalo and the connecting line between the two hosts MW and 

M31. 



Ilibeskind et al.l 120071) as we have stacked the distributions from 
two hosts. 



4 DISCUSSION AND CONCLUSIONS 

Motivated by recent claims that the two Magellanic Clouds were ei- 
ther formed close to M31 (and ejected from it by a recent merger) 
or simply accreted f r om the outer reac hes of the Local Group 
dvan den Berghll2010l : iBusha et alJuOlOf) we examined the likeli- 
hood of this situation in Constrained Local UniversE simulations 
of the Local Group. We find that subhaloes disloyal to their first 
host exist, and called them renegade subhaloes. They are subhaloes 
accreted by one host halo at a given time, yet at z = are found 
within the virial radius of completely different host. 

We have examined various physical properties of the rene- 
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gade population and contrasted them with the full loyal popula- 
tion. We find that significant differences are seen in the relative ve- 
locities, sphericity, radial distribution and spin parameter. Further- 
more, renegade subhaloes are significantly more anisotropically 
distributed than the full subhalo population - more renegades are 
found close to the line connecting our two main galaxies than not. 

Despite the differences highlighted above, detecting renegade 
subhaloes is only possible if they leave an imprint on the full sub- 
halo population - its not enough to have vastly different properties. 
One thus has to compare the dashed with the dotted line in the fig- 
ures presented above, to gauge any observable signature. In princi- 
ple, an observer will only have access to the set of visible satellite 
galaxies and hence could generate a distribution akin to the ones 
presented in Fig|4] But only if loyal satellites show a substantially 
different behavior to the combined sample will there be a chance 
to (observationally) confirm their existence (as, for instance, the 
difference in the velo city distribution for backsplash and infalling 
subhaloes reported bv lGill et al J 20051) . We thus conclude that their 
(observational) detection will be difficult unless orbital information 
can act as a discriminator. 

The abundance of renegade subhaloes appears to be a func- 
tion of both the distance between the hosts "sharing" them (which 
also corresponds to redshift in our case) and the mass of the hosts. 
Yet despite their existence, it is still a great challenge to explain 
the origin of the MCs as renegade subhaloes. In our (dark matter 
only) simulation, renegade subhaloes are not massive enough to 
realistically be called "MC"-like. Our inability to find a MC-like 
renegades, however, may simply be due to the difficulty in produc- 
ing such massive satellites in the first place. Further, our simulation 
currently only considers dark matter whereas it has been shown be- 
fore that baryonic physi cs will have a certain impact upon the prop- 
erties of subhaloes (e.g. lLibeskind et alJl2010] : lRomano-Diaz et al.l 
120080 . Future work that quantifies the frequency of finding MC-like 
objects in MW sized haloes, and extends this work by then examin- 
ing the likelihood that these subhaloes become renegades is needed 
to verify or falsify this formation scenario. 
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